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l-_ ' Abstract 

o ; 

7— ( , The extreme eigenvalues of adjacency matrices are important indica- 

T-H . tors on the influences of topological structures to collective dynamical 

kJ ' behavior of complex networks. Recent flndings on the ensemble average- 

. ,—t , ability of the extreme eigenvalue further authenticate its sensibility in the 

ji^ ■ study of network dynamics. Here we determine the ensemble average of 

$>H ' the extreme eigenvalue and characterize the deviation across the ensemble 

through the discrete form of random scale-free network. Remarkably, the 
analytical approximation derived from the discrete form shows significant 
improvement over the previous results. This has also led us to the same 
conclusion as [Phys. Rev. Lett. 98, 248701 (2007)] that deviation in the 
reduced extreme eigenvalues vanishes as the network size grows. 

Network extreme eigenvalues are succinct descriptors of the in- 
fluence of the underlying topological structure of a complex network 
on its dynamics. This makes them important predictors of epidemic 



C^ 



threshold of infectious diseases that propagate within real world com- 
plex network. Indeed, the recent demonstration that these eigenval- 
ues are ensemble averageable has provided further support for this 
view. In this paper, we study into the ensemble averageability of ex- 
treme eigenvalues through a new perspective: the connection between 
multi-modal network and scale free network. The discrete nature of 
the multi-modal network has allowed us to arrive at an improved an- 
alytical expression of the extreme eigenvalue for scale-free networks. 
The extreme eigenvalues calculated from our analytical expression are 
found to closely correspond to those obtained numerically, thus mak- 
ing significant improvement over earlier versions. The implication is 
a more accurate estimate of epidemic threshold, which is important 
for the elucidation of how vulnerable a particular network structure is 
to epidemic spreading. Our results are also applicable to the evalua- 
tion of strategies that aim to contain the spread of infectious diseases 
through the adjustment of the topology of network structures. 

Many concepts in network science have been well recognized as fundamental 
tools for exploring the dynamics of complex systems. In particular, scale-free 
networks are used widely to describe and model social, biological and economic 
systems [U HJ [31 IH [S] . In an ensemble of scale- free networks, although the de- 
gree distribution of the nodes remains the same, the topological structure of 
each individual network can be diverse with different connections introduced 
between the nodes. Such structural diversity can lead to discrepancy in dy- 
namics of the individual network. Since the structural influences on certain 
dynamical processes are governed by the extreme eigenvalues of the network 
adjacency matrices [HI El HI IB UHl HH HH US]: deviations in the extreme eigen- 
values in network ensembles are of increasing interest. Recently, it is found that 
the extreme eigenvalues of adjacency matrices, despite fluctuate widely in an 
ensemble of scale- free networks, are well characterized by the ensemble average 
after normalized by functions of the maximum degrees |14| . Specifically, it has 
been proven that the probability of having greatly deviated extreme eigenvalues 
in the ensemble diminishes as the size of the network increases. Considering the 
rich assortment of possible structural configurations for scale- free networks in an 
ensemble, this averageability is significant as it implies that dynamical processes 
which are governed by the extreme eigenvalues can be described simply using 
the ensemble average without the need of incorporating the connection details 
of the individual network. In particular, the average of network synchronization 
ability and epidemic spreading threshold are shown to be well approximated by 
functions of the ensemble average of the eigenvalues. Therefore, finding a way 
to determine the ensemble average of the extreme eigenvalues becomes crucial 
to uncover the topological influences of the network structure on a number of 
network dynamical processes. 

To the best of our knowledge, the extreme eigenvalue of adjacency matrix 
of random, undirected scale-free network has been analytically approximated 



up to the second order correction as Xj^ « kn + k)j' — 1 which gives better 
precision over the previous result Xjj « kn [El HSl [TTl [18] . Note that kn is 
the largest degree of the network and k}j denotes the average degree of the 
first nearest neighbors of node H. The probability distribution of the largest 
degree Pdikn) is given by the Frechet distribution and the ensemble average of 
kn can be calculated from Pdikn)- However, as we shall show later, both of 
these approximation to Xh can be further improved. 

In this paper, we investigate the extreme eigenvalue of undirected scale- 
free network through its discrete form, the multimodal network. For directed 
network, the extreme eigenvalue can be obtained from Rcfs. [IH |2Q]. Bene- 
fited from the mathematical properties of multimodal network that are more 
tractable, we found a way to analytically determine the ensemble average of the 
extreme eigenvalues while investigating the circumstances under which individ- 
ual network can be better represented by the ensemble average. In addition, 
for bimodal networks which are shown to be more robust than the scale-free 
networks against both random and target removal of nodes [HI [52] , we study 
the difference between them and scale-free networks in terms of the ensemble 
average of the extreme eigenvalue. 

A multimodal network |23j with m modes contains m distinct peaks in the 

degree distribution: P{k) = Y^^i i^i^{k — ki). Note that S{x) is the Dirac's delta 

function, r^ — ria^'*"^-* and ki — kib~'-'^~^^ ior i = 1,2,- ■ ■ ,m. It is assumed 

that a > 1 and < 6 < 1 such that the degree distribution of the multimodal 

network follows a power law P{ki) = r^ ex fc~ , and hence ri > r2 > ■ ■ ■ > r,„ 

for ki < k2 < ■ ■ ■ < km- As m — >■ oo, the multimodal network converges to a 

scale- free network. The largest degree of the network is km, the smallest degree 

1 
is ki which is between 1 and {k), and we have h ~ {ki/km) '""^ . The rest of the 

parameters can be determined through the following equations: 



Y^n = ri5:a-(-i) = l, (1) 

m m 

Y^hn = fciri^(a&)-(*-i) = (fc). (2) 

1=1 i=l 

We shall follow the method outlined in [TS] to find \h of the multimodal net- 
work. Let G be a graph with vertex set V{G) = {iii, W2, • • • , Wm} and adjacency 
matrix A. For each positive integer n, the number of different Vj — Vi walk of 
length n in G, denoted by yj^i{n), is the (j, i)-entry in the matrix A^. Note that 
two u — v walks W = {u — uo,ui, ■ ■ ■ ,Uk ^ v) and W' = {u = woj^ij • • • ,vi = v) 
in a graph are equal if fc = ^ and Ui = Vi for all i with < i < k |27| . For exam- 
ple, Wi — ivo,vi,V(i,V2) is different from W2 — (t^Oj^2,'yo! ^2) and both should 
be included in the calculation of yo-^2(3). In other words, yj-^i{n) is the to- 
tal number of all possible walks of length n from node j to i, including those 
go backwards. Take a fully connected network with 3 nodes as an example, 
there are 3 walks of length 3 from node 1 to node 2, i.e. 1 — > 2 — > 1 — > 2, 
1 ^- 3 -^- 1 ^ 2 and 1-^2^3-^2. This corresponds to {A^)i2 == 3. In the 



eigen-decomposition form, we have A" — vD"v' , where v is the square matrix 
whose columns are the eigenvectors of A, v' denotes the inverse of v and D is 
the diagonal matrix whose diagonal elements are the corresponding eigenvalues, 
i.e. Da — A;. Hence, 



y,^,(n) = (A"),,=VArx;,v<,. (3) 



Note that Eq.(|31) gives a summation over the nth power of the eigenvalues. When 
n is sufficiently large, the nth power of Xh will be much larger than the nth 
power of the rest of the eigenvalues. Therefore, yj^i{n) can be approximated 
in terms of only the maximum eigenvalue Xh as 

Vj^tin) ^ X'}fVj^HVH^^ . (4) 

Now if we consider the number of walks of length n+2 which start and terminate 
at node H, 

(5) 



(6) 



The first term on the right hand side of Eq. ^ corresponds to the number of 
the nearest neighbors of node H, i.e. the largest degree of the network kn- In 
|15| . the second term on the right hand side of Eq. (jB)) is shown numerically 
to be very small for scale-free networks and is hence neglected. Since we are 
interested in finding a better approximation to the ensemble average of the 
maximum eigenvalues, we shall include the second term on the right hand side 
of Eq. © in the calculation of Xh through a statistical approach. 

Starting from node H, the total number of all possible walks of length n to 
any node in the network is 

VHin) = Y. y"^^(''^ = E(^")^^ ■ (7) 

Since out of the N{k) total number of in-links and out-links, kj are directing into 
node j, hence, the fraction of these walks that end up at node j is approximately 
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Therefore, 



N{k) 



and 



k 
VH^j {n) w yn (n) -j^ , (8) 

yn^H (n) ~ yn (n) j^ . (9) 



From node j, the number of possible one-step walk is equal to the number of 
neighbors of node j, i.e. yj{l) = kj. Similarly, from a neighbor of node j, says 
j'l, the number of one-step walk is equal to the number of neighbors of node j'l, 
i.e. j/ji(l) = %i . Walking two steps from node j is the same as walking one 
step from the neighbor of node j, hence, 

%(2) = E%,(1) 

9=1 
9=1 

= fc,fcf\ (10) 



„(1) _ ^fc, 



where fc ■ ' = J2q=i ^jq/^j^ ^^ ^^^ average degree of the first nearest neighbors of 
node j. Since node j is one of the neighbors of nodes jq, thus among the two- 
step walks that begin from node j, all walks that go from node j, through its 
neighbors, and then go back to node j, is included. This means that backward 
walks are included in the calculation of yj{2) in Eq ((TO)) . Therefore, 

,,-.^(2)«fc,fcW^. (11) 

Since kjk, can be small, the approximation in Eq. (jlip may not be precise 
for each j. Hence, the approximation is applicable only as an average over all 
nodes with degree kj instead of each individual case. 

For multimodal scale-free network, there is a finite number, m of distinct 
degrees ki, each with probability r^. Thus, 



(1) 

xl^k^+y2R^kf^, (12) 



where 



I -I - 1/iV for z = m . ^ ' 



Equation p2|) implies that \h depends on the specific way the nodes within the 
network are connected, which can deviate broadly across the ensemble. When 
the exponent, /3 of a scale- free network is small, we have a more heavy-tailed 
degree distribution. This results in a larger variation in the distribution of k\ 
in the network ensemble. Hence, the values of \h in the ensemble deviate more. 
In general, k^ (x y^{k)N. For multimodal network with fixed fci, the parameter 
b and hence ki will be fixed. In order to have larger (fc), the fraction of large- 
degree node has to be higher and the fraction of small-degree node has to be 
lower, this results in a more heavy-tailed distribution with the exponent /3 to 
be smaller. In other words, /3 oc l/(fc). When the network size is larger, (k) 
has to be smaller for a fixed value of km , /3 is hence larger and the probability 



of having larger-degree node drops rapidly. Therefore, the values of Xh in an 
ensemble of multimodal networks deviate less as the networks become more 
sparse. In an ensemble of sparse networks, the individual network can thus be 
well represented by the ensemble average. 

On the other hand, for a fixed value of (fc), the degree distribution of mul- 
timodal networks vary with different values of ki. Specifically, /3 ex fci. For 
two multimodal networks A and B of same size but different values of fci , the 
fraction of large-degree node for the network with smaller ki , says network A, 
has to be larger in order for it to have the same average degree with network B. 
Hence, Pa < Pb- The choice of different values of ki can thus lead to deviation 
in Xh- Specifically, ki = 1 gives the extreme eigenvalue that is the largest, and 
Xh decreases as fci increases (see Fig. 1). 

For a random network, the average of the sum of the nearest neighbor de- 
gree is Z2 = G'q{1)G[{1) HH]. Note that Goix) = J^'^^oPk^'' is the gener- 
ating function for the probability distribution of node degree while Gi(x) = 
J^kLo kpkx'^ /{k) is the generating function for the distribution of the degree of 
the vertices which we arrive at by following a randomly chosen edge. Hence, 
for a random multimodal network without any degree-degree correlation, k^ ' « 
{k^)/{k) and 
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Note that the second moment (fc^) is a converging function of to, specifically, 

with (afe^)^™ — > for TO — > oo. 

With Eq. [T31 we study the dependence of the ensemble average of the 
extreme eigenvalues on the mode number to, for multimodal networks. Here, 
we set fci = {k) /2 and km — y^{k)N. Figure 2 show the ensemble average 
of the extreme eigenvalues for multimodal network with to. modes. {Xh) is 
the highest for the bimodal network and it decreases gradually as to, increases 
and converges to a finite value. Meanwhile, it can be seen from Fig. 2 that 
(Xh) is larger for bimodal networks. Thus, although the tolerance against both 
random and targeted removal of node is optimal for bimodal network, epidemic 
spreading is less controllable. 




Figure 1: Dependence of (Xh) on fci for the multimodal network with (fc) — 6, 
AT = 3 X 10^ and m = 2 (solid line), 10 (dashed line) or 21 (dotted line). 
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Figure 2: Dependence of (Xh) on m for the multimodal networks with fci = 
(fc)/2. Note that sizes and average degrees of the networks are: (1) (k) = 4, N = 
3000 (circles), (2) (fc) = 6, A^ = 3000 (squares) and (3) (fc) == 6, A^ = 8000 
(triangles). 
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Figure 3: Dependence of (\h) on (a) (fc) and (b) N for ensembles of scale-free 
networks with fci = (A:)/2. Note that A^ = 3 x 10 ^ for (a) and (fc) = 6 for (b). 
(Aij) = \fkH are shown in solid curves, (A^) = \JkH + (fc^)/(fc) — 1 are shown 
in dotted curves, analytical results from Eq. (1141) are shown in dashed curves 
and numerical results of the B A model averaged over 200 realizations of network 
are shown as squares. 



To verify the accuracy of Eq. [U, we compare analytical results of Eq. [HI to 
numerical results. For this, we generated a scale-free network with k\ = {k)/2 
and km = \/{k)N using the Barabasi- Albert (BA) model [2]. An ensemble 
with randomized network topology is then created using the degree-preserving 
algorithm of Ref. [29]. For uncorrelated networks, we choose networks with 
assortativity coefficients near to zero. The maximum eigenvalue of each network 
is then computed and an ensemble average is obtained. In Figs. 3, we show the 
dependence of (Xh) on (fc) and N. Note that the numerical results are shown as 
squares. Next, we compute the ensemble average for multimodal networks with 
the same parameters using Eq. (fT4|) by having r™ > 1/A^ so that there is at 
least one node with degree km- In addition, we compare our results with those 
/kn and (Xh) 



given by {Xh) = V^h and {Xh) = \JkH + (fc-^)/(fc) — 1. As shown in Fig. 3, 
our results give values of Xh that are closer to the numerical results compared 
to approximation from the previous results. 

Many real- world networks are not uncorrelated, instead they show either 
assortative or disassortative mixing on their degree. For instance, the physics 
coauthorship network in Ref. |30) is assortative while the world-wide web net- 
work is disassortative [2] . For ensembles of network with identical degree distri- 
bution, {Xh) of assortative networks with a preference of high-degree nodes to 
link to other high-degree nodes are larger than {Xh) of disassortative networks. 
For these networks, k\ ex k~^ with v > ^ for disassortative networks and ^ < 
for assortative network. Hence, although Eq. ([Tl)) gives ensemble average of Xh 
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Figure 4: Numerical results for (a) the distribution P(A^) of the normalized 
extreme eigenvalue A^ for N = 1000 (dotted line), 3000 (dashed line), 4000 
(solid line) and (b) the N dependence of the corresponding standard deviation 
a^ . Note that all ensembles of network consist of 5000 realizations of the 
networks with (fc) = 6 and ki = 3. 



for randomly connected networks, it can be generalized as 
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for correlated networks. In fact, deviation in the extreme eigenvalue is larger in 
network ensemble with varying assortivities. Nonetheless, as shown in Ref. [14] . 
fluctuation in the normalized extreme eigenvalue diminishes as the network size 
grows. In Fig. 4, we show the distribution of the normalized extreme eigenvalue 
A^ for iV = 1000, 3000 and 4000. Note that a BA network is first generated, then 
an ensemble is generated by implementing (J^i^i)"^ links rewirings following 
the constraints outlined in Ref. [29] and degree correlations in the networks 
are those generated due to these constraints [14]. As the network size grow, 
the distribution of A^ becomes more peaked and the standard deviation a^ 
decreases. 



In conclusion, the ensemble averages of the extreme eigenvalues of scale-free 
networks can be determined more precisely through the multimodal networks 
with a large number of modes. Previous approximations give much lower values 
on ensemble average of the extreme eigenvalues, and this can cause an over- 
estimation of epidemic threshold. When dealing with network dynamics such as 
the epidemic spreading of the community-acquired meticilin-resistant Staphy- 
lococcus aureus (CA-MRSA) superbugs that are resistant to many antibiotics 
[31] . over-estimating the epidemic threshold can lead to serious consequences. In 



view of this, the analytical solution derived from the multimodal network which 
is able to provide significantly closer approximation to the ensemble average of 
extreme eigenvalue of scale-free network is important. We have demonstrated 
that our analytical approximation predicted accurately the ensemble average of 
the extreme eigenvalues for scale free networks with /3 « 3 and km = ^jN{k). 
In fact, Eq. ([H]) is valid for a broad class of scale- free networks with different 
values of /? and km- While km is a free parameter, the exponent j3 can be ad- 
justed by tuning parameters a and h through the relation: /? = — In a/ In 6 -I- 1 
[23] . From Eq. (fT4|) . it is clear that {Xh) increases with an increase of \/km- 
Furthermore, it can be deduced from Eq. P^ that {Xh) decreases as /3 in- 
creases. This results from a drop in the variance in degree and the fraction of 
high degree nodes, as the exponent /3 increases. 
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